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I. INTRODUCTION 



The physics related to the IJ=00 vrvr and KK coupled-channel system, especially the 
property of the narrow /o(980) resonance, has been a subject of long and wide interests 
both experimentally and theoretically (For an incomplete list of references, see for example, 
Ref. [0 - |2^). In the latest version of Review of Particle Physics ||2^, the width of the 
/o(980) resonance is estimated to be about 40 to 100 MeV and it is remarked that the "width 
determination very model dependent". Previous estimates on the width of the /o(980) 
resonance may essentially be categorized into two classes: the first is the S matrix or T 
matrix pole fit, and the second is from the invariant mass fit in production experiments. 
It should be stressed that in a production experiment the spectrum is highly sensitive to 
the structure of the production vertex and therefore the invariant mass fit should not be 
considered as very accurate. A numerical calculation to the formfactor-like quantity in the 



coupled-channel system reveals that the behavior of the coupled-channel form-factor 
around the /o(980) region is very fragile: it can easily generate a peak or a dip, and the 
difference only comes from the (bare) production vertices. For the 5* matrix fit or the 
T matrix fit, in the literature various models have been proposed to determine the pole 
position of the /o(980) narrow resonance. Different from these approaches we in this paper 
will use the dispersion techniques to study the IJ=00 vrvr and KK coupled-channel system. 
The theory of dispersion relation starts directly from fundamental principles like unitarity 
and analyticity and when combining with experimental information the dispersion theory is 
predictive. As we will see from the discussion in this paper, the analytic structure around the 
/o(980) or the KK threshold region is very complicated: there are branch point singularities 
which are very close to the /o(980) pole, and it is out of question that one has to carefully 
examine the influence of these nearby branch point singularities to the determination of the 
pole location of the /o(980) resonance. Furthermore, the vrvr scattering phase shift 5^, as 
an analytic function defined in the single channel unitarity region, is different from the 5^ 
defined in the coupled-channel unitarity region and the two quantities have to be treated 
separately in using dispersion relations. Though the analytic structure around the /o(980) 
region is very complicated, physics related to the /o(980) resonance can be discussed by 
using properly established dispersion relations. 

For the simpler case of vrvr scatterings in the single channel approximation, we have in 



Ref. set up a dispersion relation for the real part (defined in the physical region) of the 



T matrix which is related to the experimental observable sin(257r) by the following equation, 

in which F = T(l + S)/S and F is equal to 2ReT in the physical region. In Eq. (|lD the 
position of zeros of 5* on the physical sheet of the complex s plane are denoted as Zi. The 
discontinuity of F on the left-hand cut {l.h.c), L = (— oo,0] manifests itself in the left- 
hand integral on the r.h.s. of the above equation, where one subtraction to the integral is 
understood. A very attractive feature of Eq. (|I]) is that it explicitly shows the contributions 
from different types of dynamical singularities: the resonances and the left-hand cut (and 
bound states or virtual states can be easily included when necessary). Eq. (|l]) has proved 
to be useful in determining the pole position of the a resonance p9i , especially in clarifying 
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the role of the left-hand cut in the determination of the a resonance. That is, the l.h.c. 
contribution to s\.n{25.,^) is negative and concave and therefore there must exist a a resonance 
to explain the experimental data. 

The aim of this paper, as stated above, is to extend the discussion of Ref. p9| to the more 
realistic case of the IJ = 00 vrvr and KK coupled-channel system. We will set up dispersion 
representations for physical observables in which the unitarity cut contributions are dissolved 
and the dependence of a physical quantity on the left-hand cut [l.h.c) becomes explicit. The 
method proposed in this paper is rather general and can in principle be extended to more 
complicated situation with more than two channels as well. 

The material of this paper is organized as follows: The Sec. |I| contains the theoretical 
discussions on the coupled-channel system. In Sec. [II A| the general property of the analytic 
structure of the tttt and KK coupled-channel system is reviewed, including the analytic 
continuation and the left-hand cuts. In Sec. |II B| , dispersion relations are set up for analytic 
functions constructed from the scattering T matrix. Such analytic functions contain no 
unitarity cut. In Sec. [11 C| , the analytic properties of the newly constructed functions are 
analyzed in perturbation theory and in a simple narrow resonance model. The Sec. |TT| is 
devoted to physical discussions using the dispersion relations founded in Sec. In Sec. pi A| 
the physical interpretation of the dispersion relation, Eq. (|TB[), above the KK threshold is 
discussed. It is shown that the tttt scattering phase 6tj, defined above the KK threshold 
contains a branch point singularity at s = Amj^ — 4m^. In Sec. |IIIB| the dispersion relation 
for sin(257r) in the single channel unitarity region is established which is a generalization of 
Eq. (|l|) in the coupled-channel situation. The Sec. |^ is for the application of our method in 
phenomenology where the property of the /o(980) resonance and the influence of the nearby 
branch point singularities to the determination of the pole location of /o(980) is carefully 
examined. In Sec. [lV A] a general discussion is given to /o(980) and its nearby branch point 
singularities and Sec. [IV B| is devoted to fitting the pole position of the /o(980) resonance 
combining the single channel and the coupled-channel dispersion relations set up in this 
paper. In Sec. |IV Q a subtlety in our procedure is discussed. The Sec. is for the final 
conclusion. 



II. THE ANALYTIC STRUCTURE OF vrvr AND KK COUPLED CHANNEL 

SYSTEM: GENERAL DISCUSSIONS 

A. The analytic continuation on the Riemann surface 

As is well known, the coupled-channel partial-wave scattering matrix, S, is discontinuous 
on the unitarity cut. In the two channel approximation the unitarity cut includes two 
kinematical singularities: the first one starts from 4m^ to oo whereas the second cut starts 
from Arrij^ to oo, which defines a four-sheets Riemann surface. The relation between the S 
matrix and the T matrix is, 

^11 = 1 + 2tpiT(, , 
522 = 1 + '2ip2T22 , 
5*12 = 2z^piP2T/2 . 

(2) 
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For the T matrix the unitarity relation on the unitarity cut reads, 



ImTn = TnpiT*! x d{s - Ami) + x e{s - Ami) , 

ImTu = TuP2T;, X e{s - Am\) + TnPiT*2 x Q{s - Am^^ , 
Imr22 = T22P2T2*2 X Q{s - Am\) + T^-,p{ri^ x Q{s - Ami) , (3) 



where T12 = T21 has been used. The kinematic factors pi and p2 are, pi = Jl — Am^/s 



and p2 = yl — 4m|^/s, respectively. One attempts to extend Eq. (j^) down to the lowest 
threshold. That is, below the second threshold, the above equations may be written as, 

ImTn = TiipiTl^ , 
Imri2 = TiipiTi*2 , 

ImT22 = T21P1T12 . (4) 

But in general Eq. (Q) would not be correct in the presence of the anomalous threshold and 
when the left-hand cut comes across the Ami threshold along the real axis. In the present 
case of TTTT and KK system there is no anomalous threshold, but the latter does happen. 
In the KK KK amplitude the left-hand cut starts from 4m|f — Ami to —00. As a 
consequence the third equation of Eq. (|[) only holds true above Am\ — Ami along the real 
axis rather than above Ami. The l.h.c. in Tn and T12 is (—00, 0], and L = {—00, Am\ — Aml] 



for T22 (for more information on the location of the l.h.c, see Ref. ||30|| ). 

From Eq. (^) the analytic continuation of T or to the second Riemann sheet can be 
made, Q 

rpll _ 1 rpl 

1 + 2iT{^pi 

rpll _ 1 rpl 

1 + 2zT/iPi 

^22 = ^22 — 2iT2iPlT(2 , n-rpl ) (5) 

1 + 2tTlipi 




= T^Bii = l+2ipiTl, l+2ip,Tl, . (6) 



or in short notations. 



From Eq. one also obtains the T matrix on the third and sequentially on the fourth 
sheet. 



rpUl _ rrilpi^^^ _ rpl _ rpl / detS detS 1 ^'7^ 

= ™ ~ 1 + 2zTlp ~ I -^^P^^i l+2ip,Tl, I ' KV 

^ ^ detS detS ' 



^We simply point out that, assuming real analyticity, the space between 4m — 4m^ and 4m^ 
already allows one to make analytic continuation to the second sheet. One may make use of the 
Kallen and Wightman theorem to refine the analysis, but in here we will not discuss the subtlety, 
the content goes beyond the scope of this paper. 
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T^v ^ X^B/v = TM ^2^1^ ? ) . (8) 

Similar discussions on the analytic continuation of the hadron form-factor, F = {Fi, F2), 
can also be made. The difference between F and T is that the former is free from left-hand 
singularities on the physical sheet by construction. The spectral function of F is obtained 
through the LSZ reduction formalism, 

ImFi = FipiT*! + F2P2T*2 , 

ImFs = ^2^2^2*2 + F.piT*, . (9) 

The analytic continuation of F is similar to Eqs. (il)-(i), 

F^^ = F^Bji , F^^^ = F^Bjii , F^^ = F^Bjy . (10) 

The continuation to the fourth sheet can either be obtained by analytic continuation from 
the second sheet or from the third sheet, which gives of course the same result. 



B. The dispersion representation with only left hand singularities 

To proceed we first notice that the method used in Ref. [^] to derive Eq. (|ip, which 
requires analyticity of the spectral representation, is not applicable in the coupled-channel 
situation. For example, one may set up the coupled-channel dispersion integral equation for 
the form-factor F on the physical sheet (throughout the text dispersion integrals are always 
written in the un-subtracted form, but possible subtractions are understood), 

TT J4m2 S' — S 

which defines a 2-dimensional singular integral equations with Cauchy kernel, where the 
integral kernel K, from the above discussion, can be written as, 

K(.) = (^^^'' P'l'^)i9is-Aml)-eis-Aml))+^P^'^jj ^^^C ) 0(. - 4m|) . 

(12) 

Apparently, K is not an analytic function on the real axis even though T is, hence the 
method in Ref. cannot be used here. Numerical solutions of Eq. (|TTp are searched for 



in the literature [|28[, but the numerical solution is not unique |^ as a reflection of the 
general mathematical theory [^, and the attempt in picking up the fundamental solution 
from others becomes extremely difficult numerically. 

However, as will be shown in the following, the analytic structure of F defined in Eq. (^ 
can be studied even though we do not know how to solve Eq. (^) analytically. The crucial 
observation is that, according to the way of analytic continuation, we have the following 
identity, 

±fi£!(ll + £ll(ll + Elll(ll + E!X(llyir^O. (13) 

2ni Jc T — Z T — Z T — Z T — z 
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where the contour C encircles the right-hand cut R =(4m|- — 4m^, oo) and along the 
contour the complex cut plan is on the left (the contour C would have encircled the entire 
unitarity cut if there were no l.h.c intercrossed the right-hand cut). Since functions Fi - 
Fiv are analytic on the entire cut plane except for isolated singularities and except for the 



left-hand cuts the integral contour in Eq. (|13[) can be deformed, and F can be expressed as 



sum of poles and left-hand integrals. The form-factor may contain bound state poles on the 
physical sheet which correspond to the bound state poles of Sj, and may contain resonances 
from Bii to Biy on the other sheets (The virtual state, if exist, can also be incorporated 
into our formalism easily, but we do not discuss them here for simplicity). The analytic 
expression of the scalar form-factor on the physical sheet can therefore be obtained. 



F{z) 



\-^^'^^r(.\^ \- ^(4"^)B(n)(4"^) , 1 fDtsciFCl 

Z Sj n=IIJIIJV Z- Z) ^71^^ JL Z Z 



/C{z) 



(14) 



where L = (— oo,4m|^ — 4m^]. In Eq. (|T3|) Sj denote the positions of bound state poles of 
F on the physical sheet and z(' , z^ and z^^ denote the position of the resonance poles on 
the second, third and fourth sheet, respectively. As can be seen from Eqs. (||)- (||), they are 
also zeros of Su, detS and S22, respectively. The notations F and B in the above equation 
denote the residues of the corresponding functions at their pole positions. The denominator 
C in the above equation is, 

2 + 517 + ^ -2^Pi^/2 + ^ 

(15) 



C(z) = 1 + Bji(z) + Bjn(z) + Bjv(z) 



Eq. ([IID provides an analytic expression of the form-factor F, and the solution of the integral 
equation (|TT|) provided that the behavior of the form-factor on the left-hand cut is known. 
However, these formulae are highly non-trivial. Since Fj is free from any left-hand singu- 
larity by construction (though F on other sheets does contain the left-hand singularities) 
and contains no resonance pole on the physical sheet, the residues of F at poles are not free: 
they have to be arranged in such a way that on the r.h.s. of Eq. ([l^ , the l.h.c. and the zeros 
in the numerator (the term inside the square bracket) and in the denominator, cancel each 
other.0 Though Eq. (0) is rather complicated, one simple way to check the correctness of 
Eq. ([T^ ) is to look at the case when the coupled-channel decouples, that is T21 = T12 = 0. 
In such a case it is easy to check that the effect of the third sheet is reduced to that of sheet 
II and IV and, 

C(-) = 2 ^ . (16) 

It is then easy to figure out that Eq. ( |I^ ) is reduced to two independent single channel 



expressions 29 



^The cancellation has only been verified in some very simple cases, see Ref. 
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In principle, our method presented above can be generahzed to the case with more than 
two coupled-channels, since the key point in deriving the analytic expression, Eq. ([T3|) , can 
be extended to higher dimensional case. Also, Eq. (|T^) needs more detailed analysis both 
theoretically and phenomenologically. But we will not discuss these topics here. Instead we 
will turn to discuss the analytic property of the coupled channel T matrix. 

The above method can be easily applied to discuss the analytic structure of the T matrix 
itself, since the only difference between T and F is that T is itself discontinuous on the left. 
We have, 

1 r T^(t) T"fr) T™fr) T^^(t) 

Zm Jc T — Z T — Z T — Z T — z 

from which we get 

1 /• Disc fx^C) 
{z) C (z) + $ {z) - — i z^^^dr = , (18) 

where 



27rz Jl t — z 



(TC)n = 7;-\Su- — -i^ + ^) 

detS\ 



1 




1 


'S'22 


2ipi 


Sii 


detS 


1 


— 


1 


'S'li 


2ip2 


'S'22 


detS 



(TC)22 - 77— 'S'22 - -7-77; + 



Sn J 



(TC). = T.ll + i- + i- + -^) . (19) 



and 



Res[TC(s,)] ^="^'''^Res 



TC(z] 



In Eq. (|20|) the first term corresponds to the sum of bound state poles and the bound state 
may couple to channel 1 or 2, or both. The residues can be read off from Eq. (|T^. The 
second sum on the r.h.s. of Eq. (^) corresponds to the contribution of resonance poles. It is 
worth noticing here that not only in (TC)22 but also in (TC)^^ and (TC)-^2 the l.h.c. starts 
from 4m 1^ — 4m^. One can prove that the matrix function TC defined in Eq. (|19[) satisfies 
the following property, 

TC(pi, P2) = TC(-pi, P2) = TC(pi, -P2) = TC(-pi, -P2) . (21) 

In the case of tctt, KK scattering we assume there is no bound state pole,0 so in an 
explicit form the function $ defined in Eq. (00) can be written as. 



^See however Refs. |TT|,|T9|,|T | . 
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^ _ _ I 822(2} J j)/2ipi{zjjj) ^ detS{z}y)/2ipi {z}y) 



2^pl(4,)^(l(4,)(4, -z) ^ {detSnzhj){z}jj - z) ^ S'^^{z\y){z\y 



S'ii{zh){zh -z) ^ {detSy{z}ni){z\n ' z) Y 2zp2(4y)^22(4y)(4y - z) 



C. The analytic property of the function TC 

A question naturally arise at this moment is that what we have done by transforming 
the matrix T to TC as expressed by Eq. (0)? The answer is of course that under such 
a transformation the matrix function TC is real analytic on R. The imaginary part of 
TC in the region (4m^,4m|' — 4m^) comes solely from the dynamical effects of the l.h.c, 
as TC(pi,p2) = TC(— pi,p2) in this region. To see this more clearly one may perform a 
perturbation expansion, for example a chiral expansion, on TC. Up to O(p^) term, one may 
write, 

(TC)n = ^ {8^Pl{T^? + T^t^) + Sp^T^f + Spip^T^^Ti?) + 0(/) , (23) 

where the superscripts in the small brackets denote the order of chiral expansion in powers 
of p^. Along the real axis above the second threshold, one can use the perturbation unitarity 
relation to recast Eq. (p3D as, 

(TC)ii = 4ReTi?^+^'^ + 0(/) (24) 

which is indeed real. Below the second threshold but above the first threshold one uses the 
first equation in Eq. (Q) and obtains, 

(TC)n = 4ReTi?)+(') - 4tp2TifTi? + 0{p') , (25) 



which is again real analytic. One should not draw a conclusion from Eq. (|25| ) that (TC)ii is 
real analytic down to 4m^. A close look at Eq. (|T^) reveals that (TC)ii contains the l.h.c. 
up to 4m|^ — 4m^ which comes from sheet /// and sheet IV of the Riemann surface defined 
by the S matrix, but these cuts do not show at order 0(p'*). In fact, it is straightforward to 
check the appearance of the l.h.c. when we expand (TC)ii to 0{p^), 

(TC)n = 4Tn - 4z(piTi\ + P2T12T21) - 8(p?Ti\ + P1P2T11T12T21 + PIT12T22T21) 

+^i{2p\T^^ + 3piP2rfiTi2T2i + 2pip^TiiTi2T22T2i + Ti2T2ipl{Ti2T2ipi + 2TI2P2)) 
+0(/°) . (26) 

The l.h.c. contribution firstly shows in the third term inside the second bracket on the r.h.s. 
of the above equation. The lowest order contribution it contains is ~ 0(p^)ImT22. It is easily 
understood that this term can not be cancelled by sequential terms in the chiral expansion 
since the latter are at least of order of 0(p^)ImT22. 
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For (TC)22 we have 



(TC)22 = ^ (8^P2T#^('^ + Sp^Tif + 8pip2Ti?Ti?) + 0(/) , (27) 

which is equally well written as 

(TC)22 = 4ReT2?+(')+0(p'') (28) 

above the second threshold, and 

(TC)22 = 4Rer4'^+(^) - 4^p2Tif + 4^(Imr4') - piT(?Ti?)0(/) (29) 

when s < 4mj^. The difference between here and the 11 channel is that in here the l.h.c. 
already appears in the physical sheet defined by S. For completeness we also list result for 
(TC)i2, 

(TC)i2 = 4(t(? + Tif) - Ai{p^T[fT[f + p2T[fTg!) + 0{p') (30) 

and the r.h.s. is equal to 4ReTi^2^^*'^^ + 0{p^) when s > Am\ and equals to 4ReT^^2''^*'^'' — 
4ip2Ti^2^T2^2^ + 0(/) when - Ami <s < Am\. 

From the above discussion it is realized that (TC)ii indeed contains a left-hand branch 
point singularity at s = 4m|^ — 4m^. The perturbative expansion is just used for pedagogical 
reasons. Of course the chiral expansion is badly violated in the energy region around the 
KK threshold, and the perturbative results should not be used to argue in favor of the 
smallness of the left-hand cut contributions. 

It is also helpful to understand more about the property of the function TC by examining 
how it behaves in some simple models. For example in a simple Breit-Wigner narrow 
resonance model, one may construct the S matrix satisfying unitarity: 

Tu = — + ^ 
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2ipi M| - s - i{pigi + p25'2) 



Mr- s~ i{pigi + p25f2) 



-^22 — 772 7 ; r ; (31) 

Mr- S- l[pigi + P2^2) 

where gi, g2 are coupling constants and Mr is the bare mass parameter. In this model the 
background phase is simulated by the function (j){s) with 

^ a{s) + ipip{s) 
a{s) — ipi(3{s) 

in which a{s) and (3{s) are real polynomials. It is easy to verify that the background phase 
as defined in Eq. (|32D obeys the requirement of real analyticity. 

The T matrix defined by Eq. ( PT| ) contains both the l.h.c. (from to — oo) and the 
right-hand cuts, due to the presence of the kinematic factors. However a simple algebraic 
calculation reveals that the function TC is a matrix of real rational functions, that is, it is 
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analytic on the entire s plane with only isolated singularities. In this sense the structure of 



the S matrix defined by Eqs. (31) and (|3^) is topologically trivial. This observation, which 
follows from Eq. (^) also applies to more general case of the K matrix parametrization 
when K takes the form of rational functions. Though topologically trivial, simple models as 
discussed above can still be helpful in revealing the qualitative picture of the mr and KK 
coupled-channel system. 

III. THE ANALYTIC STRUCTURE OF vrvr AND KK COUPLED CHANNEL 

SYSTEM: THE TWO 6^S 

A. The TTTT scattering phase in coupled— channel unitarity region 

In the physical region above the second threshold, there is the well known parametrization 
of the scattering S matrix. 

By using it Eq. (^) can be recasted as 

1 r?^ + 1 
(TC)n = -^^^sin2(5^, 

Pi V 

(TC)22 = -^^sin25^ , (34) 



1 


+ 1 


pi 


V 


1 


v'^ + i 


P2 


V 


(1 




\/iPlP2) 



(TC)i2 = ( cos(5, + 5k) + - cos(5^ - 5k) 



and from Eq. ([18[) we obtain. 



^2 



1 



1 I- Disc (t^C) 

V /II. 



Disc {t^C] 




cos(5^ + 5k) + - cos(5^ - 5k) = --^^^^^ I *i2(s) - — / ^ —dz | . (35) 

V (1 — ?7 ) 2 

In above equations the function $ is a rational function with real coefficients and the left 
hand integrals are also real analytic functions due to the property of real analyticity of the 
T matrix and L = {—oo, Amj^ — 4m^]. 

It is convenient to read from Eq. (|35|) that sin(2(57r) as an analytic function on the entire 
physical sheet, after analytic continuation, contains the l.h.c. starting from Amj^ — 4m^ 
to the left of the real axis. It is easy to understand that sin(257r) and t] have this l.h.c. 
separately in such a way that the l.h.c. starting from 4m|^ — 4m^ in their combination, Su, 
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cancels. The appearance of such a l.h.c. can be clearly seen from the following expressions 
for sin(257r) and 77 which are appropriate for the analytic continuation, 

sin(25^) = piJ^2 = ^{-Su - -^r]) , (36) 

Zt Tj ^11 



from which we read off, 



and 



7] = J 811822/ detS , (37) 



^ = , V''' = , (38) 



^2' = ^2 , V'' = l/v' , (39) 



which are correct when Eq. (^) and Eq. (|^) are correct, respectively. From these properties it 
is realized that on the real axis JF2 only contains the l.h.c, and r] only contains the cut in the 
region [4m|'— 4m^, Amj^] (where 77 is a pure phase) and the l.h.c. on the left. These look good 
at first glance for setting up dispersion relations for these quantities. But rj contains many 
branch point singularities (and so does sin(257r)) on the complex s plane (corresponding to 
the positions of resonance poles), which prevent the usefulness of these dispersion relations. 

In another combination of rj and sin(257r) which appears in the first equation of Eq. (p^, 
as we demonstrated before, the right-hand cuts cancel. Eq. (|35|) is the appropriate one 
for phenomenological discussions, as similar to the single channel case |2^. For example. 



one may use the experimental data of 6-,^, 6k and r] in the region above the KK threshold 
to study the properties of the widely concerned /o(980) resonance. Unfortunately, the left- 
hand integrals in the energy region required by coupled-channel analysis are very difficult to 
estimate theoretically. This shortcoming limits the phenomenological usage of our method. 
But our approach can still be helpful in studying the coupled-channel system, as will be 
discussed in Sec. |V|. 

B. The TTTT scattering phase in single channel unitarity region 

The analytic structure of the vrvr phase shift defined in the coupled-channel unitarity 
region (denoted as 5^^ later on) is remarkably different from the analytic structure of 5^ 
defined in the single channel unitarity region (denoted as 5^^^ hereafter): 

sin(2(5(i)) = p,J^ = i(5n - l/^n) • (40) 

The function JF defined in above equation only contains the left-hand cut starting from 
s = together with the ordinary unitarity cut starting from 4m It now becomes clear 



"^The perturbation expansion fails badly here, and the matrix Pade approximation destroys the 
single channel unitarity of the tttt amplitude by giving it the l.h.c. on the physical sheet 
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that the two 5-,^ are not analytic continuation to each other though they match at s 
One can then set up a dispersion relation for JF, 



Res[HzV)] 



ImiJ^ , 1 
as H — 



ds' 



Vr JAm]^ S' 



Am\. 



(41) 



where the poles are only from the second sheet in the absence of bound state poles. Eq. ( ^T] ) 



is the extension of the single channel result [2£] in the coupled-channel case. Comparing 



with Eq. (42) of Ref. |29] only the third term is new and the other two terms are the same. 
The first integral on the r.h.s. of Eq. ( ^T]) has been discussed at length in Ref. I^O], and 
it is fortunate that the second integral on the r.h.s. of Eq. (|4l| ) can be estimated from 
experiments since 



^)cos(25(2)^ 



(42) 



in the coupled-channel region. 



IV. PHENOMENOLOGY OF THE COUPLED CHANNEL DISPERSION 

RELATIONS 

A. The /o(980) pole and the nearby singularities 

From Eq. (|41|) we see that when s approaches Am\ from the lower side the value of the 
right integral will change rapidly since s = Amj^ is its branch point. It is well known that 
the branch point will develop a cusp structure to the real part of the dispersive integral. 
Also we see in the above subsection that sin(25^^'') contains a left-hand branch point at 
4m|f — 4m^. Both of the two branch point singularities are very close to the pole position of 
the /o(980) resonance. Therefore it is necessary to carefully investigate the influence of the 
branch point singularities to the determination of the pole position. In a typical K matrix fit 
dynamical singularities rather than poles are simulated by background polynomials. However 
the cusp structure of the branch point singularity should not be well simulated by a smooth 
background polynomial, especially when the background polynomial is being used to cover 
a large region of s. It is even reasonable to imagine whether such a cusp phenomenon could 
be responsible for the sharp rise of near the KK threshold.^ However, a numerical 
estimate to the right-hand integral in Eq. (^if) indicates that the right-hand cut effect is 
very weak. In Fig. |1] we plot the right-hand integral contribution to from two fits found 
in the literature 0,0, here the integration is performed from the KK threshold to roughly 
about ^/s ~ 1.5 GeV. As we see from Fig. ^that the contribution of the right-hand integral 
is in the right direction to increase sin(25^^^) rather rapidly, but the effect is too small to 
have any non-negligible influence to the /o(980) pole. This can be clearly seen by comparing 



^It is very impressive to notice that in Ref. Q, it was carefully tested up to 40 parameters in 
the background phase to check whether the background polynomial can be responsible for such a 
sharp rise. 



12 



Fig m with Fig 0. In the latter case sin(2(57r) jumps from -1 to +1 near the KK threshold 
whereas the contribution from the right-hand integral is one order of magnitude smaller. 
When estimating the right-hand integral we make use of Eq. (^) which assumes implicitly 
two-channel unitarity. Though the two-channel unitarity no longer holds true when the 47r 
channel and the rjrj channel become important, we expect the order of magnitude estimate 
remains valid. 

The left-hand integrals appeared in the coupled-channel equations, though can not be 
estimated quantitatively, can also be proven to behave mildly when the narrow /o(980) 
resonance is located on the second sheet (see Eq. (HO) and related discussions). 



B. A combined fit of the coupled channel and the single channel equations 

The single channel dispersion relation, Eq. (^I]) contains only the second sheet pole 
whereas the coupled-channel equations, Eq. ( P5| ) contain resonance poles on all sheets, i.e., 
sheet II, III and sheet IV. In the single channel approximation we have to introduce 4 



parameters for each resonance pole, two of them are related to the pole position and the 
another two are related to the residue, or the resonance coupling to JF. Now in the coupled- 
channel case each resonance pole contains 6 parameters, two of them are for the pole position, 
two are for the resonance coupling to T (for second sheet poles only) and/or $ii (the second 
sheet pole coupling to JF and $ii are the same), and the last two are for the coupling in $22- 
The pole coupling in $12 is not independent as can be verified from Eq. (0). It is annoying 
to have so many parameters associated with a resonance pole. These parameters are in 
principle correlated to each other, but no simple relation with stringent constraint between 
these terms are found. Therefore in here we approximate these parameters as uncorrelated 
in the fit. However the property of the /g^(980) narrow resonance obtained under such an 
approximation is expected not distorted much from the real situation, since the sharp rise 
of the TTTT scattering phase near the KK threshold is so deterministic, as can be seen from 
the following discussions. 

In the TTTT and KK coupled-channel system, following the conventional wisdom, we 
assume two second sheet poles. One corresponds to the a resonance, the another corresponds 
to the /o(980) resonance, in addition we assume there exists another 3rd sheet pole which 
simulate all other pole contributions.^ The right-hand integral is estimated using that 
depicted in Fig. |I| which is however found to have only tiny influence. We estimate the 
left-hand integral in Eq. (^) in the same way as in Ref. |2^. That is we either use xPT 
to evaluate ImJ-^, but with a truncated integral interval at — (here we take A= 600MeV, 
rOOMeV, SOOMeV and IGeV), or estimate IuiTl from the [1,1] Fade solution. [] 



"A shortcoming of the present method is that we can not distinguish the third sheet pole from 
the fourth sheet pole in Eq. (|35|). 

^We are not able to estimate the left hand cut at quantitative level. In a previous version we had 
fixed IvhTl using the Pade approximant, but then we realized the problem associated with the 
Pade approximation Q|. Fortunately the pole positions of the a and the /o(980) resonances are 
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In the coupled-channel unitarity region, or in Eq. (^H), however, we do not have a 
reliable method to estimate the left-hand integrals. We therefore simulate them by a group 
of (totally 3) polynomials. In here we restrict ourselves in a modest range of energy region, 
2mK < y/s < 1.2GeV, therefore the polynomials may take simple and smooth forms, i.e., 
constants or linear polynomials which further introduce a few more parameters. In the 
absence of narrow 3rd or 4th sheet pole close to the branch point singularity at Amj^ — 4m^ 
the validity of this approximation can be justified. Furthermore, the matching condition. 



sm 



(25«)|,=W =sin(2<5(^))|,=4^2 (43) 



is used to reduce one parameter in the vrvr channel. 

In the single channel unitarity region we use the combination of the CERN-Munich data 
and the low energy data from Kg^ decay , and especially the new experimental data from 



the E865 Collaboration |^5[. We take the K meson mass to be rriK = {jnK± +mKs)/2 in our 
fit. In the coupled-channel unitarity region we use the experimental data on 5^ + 5k and 
(1 — '?7oo)/4 from Refs. P^ , p7|] , and the data on 5^ from Refs. P^P^ .f] In order to extract 



the "experimental" value of TC in the coupled-channel unitarity region we have to shift the 
data to same points of s by proper extrapolation. With such a manipulation to the data, 
we are able to make the fit, see Fig. |^ and Figs. |^-^ for a typical example of the fit where 
the left-hand integrals in the coupled-channel equations are approximated by constants. In 
the fit we find that, even for a very simple parametrization form of the left-hand integral, 
there may exist different solutions - the position of the 'third sheet pole' is found to be 
very unstable. This even occurs when the 3rd sheet pole is correlated to the 2nd sheet pole 
via a Breit-Wigner narrow resonance model like that in Eq. (|3l|) . In the uncorrelated case 
the solutions for the 3rd (or 4th) sheet pole is quite arbitrary: The mass roughly ranges 
from 0.7 ~ 1.2GeV and the width ranges from very small values to roughly about 0.6GeV. 
This uncertainty may be partly due to the unsatisfactory quality of the manipulated data in 
the coupled-channel unitarity region (see Figs. but it is amazing to notice that such 

an uncertainty for the 3rd/4th sheet pole does exist in the literature where very different 
solutions are also found (see the compilation of the 2000 edition of the Review of Particle 
Physics for more information). Irrespective of the large uncertainty that the 3rd/4th 



sheet pole has, the a pole position and the /q^ pole position are however found to be rather 
stable against changes of the parametrization form of cuts and the uncertainties related to 
the first left hand cut. The results are summarized in what follows: 

M„ = 450 ~ SlOMeV , = 480 ~ 550MeV ; 
Mfii= 981 ~ 992MeV , F.// = 32 ~ 40MeV : 
a° = 0.24 ~ 0.27 . (44) 



Very similar to the what happens in the single-channel situation |^, the location of the a 
resonance is sensitive to the choice of the scattering length used in the fit. If the xPT result 



found not sensitive to the uncertainties related to our different choices of the left hand cut. 
®Such a choice of data set merely reflects our prejudice by excluding the results of Refs. [40,41|. 
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of Og = 0.220 ± 0.005 [0 is used as a constraint in the fit we find instead the following 
results (here we take the Pade solution of the first left hand cut for example): 



M„ = 465 ~ 468MeV , = 610 ~ 650MeV ; 
Mfii= 981 ~ 983MeV , T.// = 31 ~ 33MeV ; 

Jo 'Jo ' 

~ 0.221 . (45) 

From the results list above we find that the pole position of the /o(980) resonance is very 
stable. Furthermore the inclusion of the /q^(980) pole has only modest influence to the pole 
position of the a resonance by comparing with the results from Ref . The main reason 

of the stability is due to the special data point in the upper right region in Fig. from 
the CERN-Munich data: ^/s = 0.99GeV with the value of 5^ = 234.0 ± 12.3 degrees. The 
center of mass energy of this data point already exceeds 2mK± though it is located below 
the K^K^ threshold. In our fit we set the K meson mass as niK = {mK± + rnxs)/'^ which is 
above 0.99GeV, so we count this point in the single channel unitarity region (notice that this 
data point is very close to 2mK±)- Also considering the experimental error bar for the beam 
energy, the location and the value of this data point contributes the major uncertainty in the 
determination of the location of the /o(980) resonance (mainly affects the width). Manifestly 
our approach is sensitive to the experimental data near the KK threshold region. On the 
qualitative level, however, the picture that there exists a very narrow /o(980) pole on the 
second sheet should not be altered, as is determined unambiguously by the fact that the 
value of sin(257r) at the matching point (see Eq. (^)) is much closer to +1 rather than to 
-1. 



C. The subtlety in our approach: if the /o(980) narrow resonance does not locate on 

the second sheet 

As pointed out in the above subsection that the data point closest to the KK threshold 
is crucial in determining a narrow resonance on the second sheet, named as the /d^(980). 
One has to be very cautious under such a situation when the fit is very sensitive to a single 
data point. For example, if this point were removed from the data in fitting Eq. (0), or 
in other words if the value of sin(2(57r) at the matching point (see Eq. (^)) is closer to — 1 
rather than to +1, the fit program would prefer a solution without the narrow second sheet 
pole. Instead, it would give a solution roughly at Mj-n ~ 0.86GeV and Tj^n ~ 0.16GeV and 
a very narrow third (or 4th) sheet pole, M ~ 0.96 — 0.98GeV and F ~ 5 — 35MeV. Under 
such a circumstance it is the 3rd (or 4th) sheet pole combing with the l.h.c. in the first 
equation in Eq. (pSj), rather than the narrow second sheet pole, being mainly responsible for 
the sharp rise of the vrvr phase shift above the KK threshold. Though rather academic, it 
is worth pointing out that if such a scenario would happen the narrow pole position could 
not be reliably estimated within the present approach, because the third (or 4th) sheet pole 
will have a very strong influence to the left-hand integrals in the coupled-channel dispersion 
equations. To make the point more clear let us write down the following formula. 



^Previous results on this point can be found in Refs. |4^,44| 
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li)«c(TC„) = -4plA„.J^ + A.c., (46) 

where A;/ic is defined as 

\mT22 = Pi\Tu\^ + AihcOi^ml - Ami - s) , (47) 

in the single channel unitarity region. We have no knowledge on how to estimate Aihc and 
we assume it behaves normally in the sense that its order of magnitude and the behavior 
is similar to, say, the l.h.c. appeared in the single channel equation. Then it is apparently 
shown in the above formula that a very narrow resonance on the 3rd or 4th sheet near the 
left-hand branch point at s = Amj^ — 4m^ will greatly enhance the discontinuity of TCn 
and also the cusp structure of the dispersion integral of TCn provided that the sign is 
correct {i.e., to give a sharp rise, not a sharp decrease, to sin{26^^)) . The 3rd or 4th sheet 
pole itself would also contribute in Eq. ( p5D (but not in Eq. (^)), a combined contribution 
of the cut and the pole may therefore explain the sharp rise of sin(25^^^). But of course such 
a scenario is unambiguously excluded by experiments. 



V. CONCLUSION 

In this paper we have studied the IJ=00 tttt and KK coupled-channel system by using 
dispersion relations which are set up after a careful analysis to the analytic structure of the 
coupled-channel system. The effects of various cuts on the determination of pole positions 
are discussed and estimated which are found to be mild in general. Especially the role of the 
left hand cut starting from 4m|^ — 4m^ is clarified. Especially, the cut has little influence to 
the determination of second sheet poles. 

We therefore confirm, at qualitative level, the conventional wisdom in the widely used 
K matrix fit to simulate these dynamical singularities by mild background polynomials. 
To some extent, our analysis also justifies the dynamical assumption made in dynamical 
approaches (see for example, Refs. [|^,0) to partly neglect the t channel exchange forces 
by assuming them to behave mildly. 

We conclude from the discussion in this paper that there exists a very narrow second sheet 
pole with mass about 981 ~ 992MeV and width about 32 ~ 40MeV which is responsible 
for the sharp increase of the tttt scattering phase near (mainly below) the KK threshold. 
Meanwhile the mass and the width of the a resonance is estimated to be M„ = 450 ~ 
SlOMeV, To- = 480 ~ 550MeV. However, we pointed out that the width of the narrow second 
sheet pole is sensitive to the data near the KK threshold, and therefore future experiments 
with more accurate data are called for. We also find that the possible third sheet pole which 
is related to the narrow second sheet pole in the coupled channel Breit-Wigner resonance 
model can not be reliably fixed at present stage within our scheme. The reason is simply 
because only the data above the KK threshold can contribute to the determination of the 
third sheet pole and also the left hand cut starting from 4m|^ — 4m^ pollutes here. 

Our results on the pole locations of the a and /d^ resonances are in qualitative agreement 
with the results found in the literature, though the method used in this paper is rather 
different from those commonly used. In our opinion, the IJ=00 tttt and KK coupled- 
channel system is so interesting and important and it is always meaningful to study it from 
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different angles which can bring us more understanding and insight in the related physics. 
Our method may be generalized to others channels as well. But the dynamical singularities 
appeared in the case of unequal mass scatterings are more complicated. For example, in the 
case of ttK scattering we will run into the more complicated problem of dealing with the 
circular cut. 
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FIG. 1. The contribution from the right-hand integral to . 
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FIG. 2. A typical fit to sin(2(57r) in the single channel unitarity region. 
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FIG. 3. A typical fit to TCn in the coupled-channel unitarity region. 




FIG. 4. A typical fit to TC12 in the coupled-channel unitarity region. 
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FIG. 5. A typical fit to TC22 in the coupled-channel unitarity region. 
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